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ABSTRACT 

Future high precision photometric measurements of transiting extrasolar planets promise to 
tell us much about the characteristics of these systems. We examine how atmospheric lensing and 
(projected) planet oblateness/ellipticity modify transit light curves. The large density gradients 
expected in planet atmospheres can offset the unfavorably large observer-lens to source-lens 
distance ratio, and allow the existence of caustics. Under such conditions of strong lensing, 
which we quantify with an analytic expression, starlight from all points in the planet's shadow 
is refracted into view producing a characteristic slowing down of the dimming at ingress (vice 
versa for egress). A search over several parameters, such as the limb darkening profile, the planet 
radius, the transit speed, and the transit geometry, cannot produce a nonlensed transit light 
curve that can mimic a lensed light curve. The fractional change in the diminution of starlight 
is approximately the ratio of atmospheric scale height to planet radius, expected to be 1% or 
less. The lensing signal varies strongly with wavelength — caustics are hidden at wavebands where 
absorption and scattering is strong. Planet oblateness induces an asymmetry to the transit light 
curve about the point of minimum flux, which varies with the planet orientation with respect 
to the direction of motion. The fractional asymmetry is at the level of 0.5% for a projected 
oblateness of 10%, independent of whether or not lensing is important. For favorable ratios of 
planet radius to stellar radius (i.e. gas giant planets), the above effects are potentially observable 
with future space-based missions. Such measurements could constrain the planet shape, and its 
atmospheric scale height, density and refractive coefficient, providing information on its rotation, 
temperature and composition. For HD 209458b, the only currently known transiting extrasolar 
planet, caustics are absent because of the very small lens-source separation (and a large scale 
height caused by a high temperature from the small separation). Its oblateness is also expected 
to be small because of the tidal locking of its rotation to orbital motion. Finally, we provide 
estimates of other variations to transit light curves that could be of comparable importance — 
including rings, satellites, stellar oscillations, star spots, and weather. 



Subject headings: planetary systems — planets: atmospheres — microlensing — stars: atmo- 
spheres 



1. Introduction 

The exciting discovery of the transit of HD 209458 (Charbonneau et al. 2000a, Henry et al. 2000, Jha et 
al. 2000) allowed the first direct measurement of the physical properties of an extrasolar planet, including its 
radius, and together with radial velocity measurements, its absolute mass and average density. Photometric 
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followup of radial velocity planet candidates, and future ground- and space-based surveys are expected to 
discover hundreds more transiting planets. Planned or proposed space-based missions, in particular (MOST, 
MONS, COROT, the Eddington Telescope, the Kepler Mission), are expected to be able to detect flux 
variations at the 10~ 5 level. Recent Hubble Space Telescope observations of the HD 209458b transit (Brown 
et al. 2001), with an accuracy of 10 -4 , illustrate the capability of high precision space based photometry. 
With the anticipated large number of transiting planets, and the high accuracy with which they can be 
monitored, it is important to explore small variations in the transit light curve that might be detectable and 
allow us to deduce further properties of the transiting planets. 

In this paper we focus on how atmospheric refraction, which we also refer to as atmospheric lensing 
or lensing for short, modifies the transit light curve. We treat the general case of an ellipsoidal planet, 
and the same calculation also allows us to examine the effects of (projected) planet oblateness/ellipticity, 
regardless of whether or not lensing is important. We examine under what conditions atmospheric lensing 
and oblateness signatures might be significant, detectable, and distinguishable from other effects. 

We develop the lensing formalism in §^|, and contrast atmospheric lensing with gravitational lensing. For 
readers with a background in the latter, we point out several places where intuition gained from gravitational 
lensing fails for atmospheric lensing. In §^ we give a complete list of all the parameters employed in our 
model. The condition under which atmospheric lensing produces strong lensing (i.e. existence of caustics) 
is taken up in §^|. An analytic expression describing this condition is given. We emphasize the connections 
between atmospheric lensing by extrasolar planets and atmospheric lensing by solar system bodies during 
their occultations of distant stars; the latter has long been observed (see e.g. Elliot & Olkin 1996; Hubbard 
1997, and references therein). The difference is: in the extrasolar case, the lens-source distance is much 
smaller than the lens-observer distance, and the source is extended compared to the lens, while the opposite 



is true in the solar system case. In §5.1, we examine how the lensing modifications to the light curve vary 



with the model parameters. In §5.2, we address the question: to what extent can a lensed light curve be 



confused with a nonlensed light curve with different parameters? Similarly, in § |5.3| we isolate characteristic 
modifications by the planet's oblateness to the transit light curve. The issue of how absorption affects 
the color dependence of lensing is taken up in §0, where we also discuss how lensing impacts the planet 
transmission spectrum. Finally, we conclude in §ffl with a summary and with a list of other effects that might 
be of comparable magnitude to the lensing and oblateness signatures. 

An interesting paper by Hubbard et al. (2001) recently appeared as our paper was nearing completion. 
They discussed some of the same issues addressed here, but focused mainly on detailed predictions for the 
case of HD 209458b. They also presented calculations of atmospheric glow from Rayleigh scattering which 
we do not treat in this paper, although we do discuss extinction by Rayleigh scattering. 



2. The Lensing Model 
2.1. Formalism 

Much of the formalism presented here closely parallels that used for gravitational lensing, and the 



similarities and differences are discussed in §2.2. Equivalent techniques have been applied to lensing of 
distant stars by solar system bodies (Elliot & Olkin 1996; see also Draine 1998 for a treatment of gaseous 
lensing by large spherical clouds). Ray optics provides the following starting point for lensing: 
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where Op is the angle of deflection, 9s and Bj are the source and image positions, Dls denotes the distance 
between the lens and the source, and -Dol is the separation between the lens and the observer. The angle 
of deflection is determined by spatial gradient of the refractive index n: 

e D = - Vn x A£ x i, (2) 

where the sum is over individual segments of the ray A£ with £ as the unit vector, and where V is the 
spatial gradient. 

For atmospheric lensing, n = 1 + ap where p is the gas density and a is a refractive coefficient that 
depends on both gas composition and wavelength. A common parametrization, known as Cauchy's formula, 
gives ap = A\(l + Bi/X 2 ) where A is the wavelength, and A\, B\ for common gases are given in Table [j] 
(Born & Wolf 1999 p. 101). For a H 2 - He mixture with 24% He by mass, a = 1.243 cm 3 g _1 at 4400 A, and 
a = 1.214 cm^" 1 at 6700 A. 

We are interested in cases where the deflection angle is small, and the lens is thin (i.e. distances over 
which significant deflection occurs is small compared to Z?ls an d -Dol)- Suppose the axis connecting the 
observer and (the center of) the lens points in the z direction. We have 

/OO />OC 
ndz = Vj_ / apdz, (3) 
-oo J — oo 

where Vj_ is the spatial gradient in the x and y directions. We use b = (61,62) to denote the impact 
vector — the vector in x-y plane from the lens center to the point of impact. 

The inverse of the magnification matrix is 

A- 1 ee M (4) 

= 6ii + Doi + Duz8b j db i L 00 apdz > 

where we have used b = Dol^i- 

We are interested in a density profile p which is ellipsoidal in general, in the sense that p = p(r) where 
r 2 = x' 2 /a\ + y' 2 /a\ + z' 2 /a 2 . f\ The axes denoted by x', y', and z' are generally not lined up with the x, 
y, and z axes defined before. We show in Appendix A that 



• such a profile is well-motivated; 

• in the simple case of an isothermal atmosphere in hydrostatic equilibrium, the density profile is p(r) = 
poexp[— (r — Rq)/H] where Rq is a reference radius whose choice will be fixed below, po is the density 
at Rq, and H = (fc&T/ gfimn) is the atmospheric scale height where kb is Boltzmann's constant, T is 
the temperature, g the surface gravity, p the mean molecular weight, and mn the mass of the hydrogen 
atom; 

• The quantity r, which we loosely refer to as "radius" , can be written as r = yj (1 — e)b 2 + (1 + e)b 2 + z 2 
by a suitable rotation of axes, along with rescaling of the density field — as long as the deviation from 



A spherical planet is simply a special case within the class of models we study here. 
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spherical symmetry is small: e <C 1, or Oi, ai and 03 are not too different from each other. The 
projected oblateness e is related to the actual oblateness of the planet, 6a = (di — a-a) /d\ (assuming 
it is axially symmetric with ai = a 2 ), by e — e^(l — cos 2 /3) where (3 is an Euler rotation angle.0 j3 is 
the angle between the axis of rotation of the planet and the line of sight (see Appendix A for details). 
Assuming j3 is randomly distributed implies e = on average. 

We also work under the simplification that a, the density- independent refractive coefficient, is indepen- 
dent of position on the planet. This is a simplification because the atmospheric composition — hence the net 
value of a — is expected to vary with atmospheric depth. 

Putting the above together, and using Vj_p(r) = (1 =F e)(h/r)dp/dr (upper/lower sign for 61/62) in 
equation (^]), the lensing equation (equation (Q)) can be written as 

9 S = 6\ + (1 - e)0)^(u), 6% = 9j + (1 + e)^(u) (5) 

where ip(u) = a— — — / 

Dol + D LS J_oo or r 

P = Po exp[-(r - Ro)/H] 

r 2 = Dl L u 2 + z 2 , u 2 ee (1 - e)6} 2 + (1 + e)ef . 

With the above form, the problem of predicting image position(s) given a source position can be reduced to 
solving a simple single variable equation — this and associated computational tricks are discussed in Appendix 
B. 

The magnification is given by (equation ([|)): 

A ee detA i:j = [1 + 2tp + (1 - e 2 )i{? + (1 - e 2 )u 2 ^ + (u 2 + ev)^ 1 (6) 



where v = -(1 - e)6\ + (1 + e)6f 
u 2 ee (1 - e)9} 2 + (1 + e)6f 



= = aD Oh- 



udu ^Dql + D ls 



1 d 2 p 1 dp 

r 2 dr 2 r 3 dr 



dz. 



and where u, ip and p are as described in equation (g). The caustic is defined by source positions where A 
diverges. The critical curve is the image of the caustic. We discuss in Appendix B how to find both, if they 
exist — a situation referred to as strong lensing. 

Finally, the observed flux from a star during an extrasolar planet transit is given by: 

F(t) = J d 2 9 s AI(6 8 - e.{t))W(6i) (7) 

where I (9s — 9*{t)) is the surface brightness of the star as a function of source position. The symbol 0*(t) 
denotes the position of the star's center as a function of time. We have chosen the origin of 6s and 6j to be 
centered at the planet. The sum is over all images for a given source position. The kernel W(6j) describes 
the occultation and absorption; in this paper, we focus on a simple model in which W(0z) is a step function, 



W(0i) = 1 if ^ {I - e)0) 2 + (I + e)9f > R a /D OL (8) 
= otherwise. 



2 In this paper, we use the terms oblateness and ellipticity interchangeably. 
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This step function specifies Ro. ft is the (elliptic) radius below which the projected density exceeds some 
value such that the atmosphere becomes completely opaque, or it is the radius at the rocky surface of 
a planet. A more realistic treatment of absorption will have W changing more gradually than this step 
function, and also changing with wavelength; this will be discussed in §^|. As we will see, our step function 
model is actually a good approximation to reality. Note that both / and A above are functions of wavelength. 
Finally, for most purposes, we are only interested in the normalized F(t) i.e. F(t) (equation^)) divided by 
its asymptotic value well away from the transit — the stellar flux: J d 2 9sl(9s)- From now on we use F(t) 
to refer to the normalized value, and whenever we refer to "flux" in this paper, we always mean the stellar 
flux normalized by its pre- or post-transit value. 

To summarize, equations (||), (^), and (R) completely specify the problem of atmospheric lensing during a 
planetary transit, for an isothermal atmosphere in hydrostatic equilibrium with an ellipsoidal density profile. 



2.2. Atmospheric versus Gravitational Lensing 

Gravitational lensing can be described by essentially the same equations presented above, except that 
the refractive index n in equation (^) is equal to 1 — 2</>, where tf> is the gravitational potential, instead of 
1 + ap. That a is wavelength dependent implies atmospheric lensing is color dependent whereas gravitational 
lensing is achromatic. 

Similarly, gravitational lensing by an elliptic potential obeys equations (|^), (JsJ) , and (Q) with ap re- 
placed by —2cj>. It is interesting to note, however, that a well known theorem in gravitational lensing, the 
magnification theorem, does not hold in atmospheric lensing. The magnification theorem states that for a 
given source position, the magnification of all images must sum to at least unity. This can be traced back 
to the fact that V 2 </> is proportional to the mass density which is positive definite. For atmospheric density, 
the relevant quantity, V 2 p, is not guaranteed to be positive definite. In other words, in atmospheric lensing 
the magnification need not sum to unity, and thus a net suppression of flux can occur. 

Another important difference between gravitational and atmospheric lensing is that gravitational lensing 
is almost never significant when the lens is very close to the source or the observer. This arises from the 
fact that the combination of distances -Dls£*ol/(-Dol + -Dls) becomes small (dominated by the smaller of 
the two distances) if Z?ls *C -Dol or Dql <C -Dls- While the same factor applies to both types of lensing, 
atmospheric lensing has the advantage of having available an exponential density field (or nearly so: see 
Appendix A). The analog for gravitational lensing, an exponential potential, almost never occurs in nature; 
power-law fall off is far more prevalent. An exponential density profile allows large gradients to offset an 
unfavorable combination of distances. To be more precise, while a power-law profile would result in factors 
of -Dls-Dol/[(-Dol + D LS )i? ] in the relevant expressions (e.g. ^p(u) in equation (||)), an exponential profile 
gives -Dls-Dol/[(-Dol + D LS )_ff], which is considerably larger (because atmospheric scale height H <C planet 
radius Rq). Nature offers a nice existence proof: atmospheric lensing of distant stars by the solar system 



planets have been observed, in spite of the fact that Dol <C -Dls (see §0). The extrasolar case we are 
interested in can be viewed as the symmetrical opposite with Z?ls -Dql- 
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3. Parameter Accounting 

Several parameters enter into the problem of planetary atmospheric refraction during a transit. However, 
most of them appear in a few combinations. It is helpful to list them explicitly. 

The quantities ip(u) and u 2 tp(u) in the lens mapping equation (|5|) can be well approximated by the 
following expressions if uDo^/H 3> 1. This is true in realistic cases because m-Dol is constrained to be 
larger than Rq in order for an image not to be obscured (occulted or absorbed) by the planet. 

i/j(u) = -B^J^IuDol/H)- 1 / 2 exp[-( Uj D OL - R )/H] (9) 
u 2 i>(u) = B^/7/2(uD OL /H)^ 2 [l + (uDol/H)- 1 } exp[-(uD OL - R Q )/H] 

B^2a^ DhsD ° L . 
H D OL + D^s 

Note that refraction is important (or the magnification A is significantly different from unity) only if 
uDol is close to Rq. Therefore, as far as the gross lensing behavior (or lack thereof) is concerned, only three 
parameters are important, B as defined above, Ro/H, and e. 

In addition to the three lensing parameters, we have three transit parameters. They are the transit 
impact parameter (as opposed to the lensing impact parameter) #tr.imp. (the distance of closest approach 
between the center of the planet and the center of the star in angular units), the transit-impact-angle 'ytr.imp. 
(the angle between the major axis of the projected planetary ellipse and the direction of transit motion), 
and the transit velocity w in angular units per unit time. See Figure |^ for an illustration. Since w basically 
rescales the time axis in our prediction of the transit light curve, we do not treat it as a free parameter in 
our predictions. We estimate it by using w = -J GMq / Z?ls /-Dol- 

Finally, there are the parameters that describe the star: the stellar radius i?*, and the surface brightness 
of the star, or its limb darkening function, which is often parametrized as: 

I{0 a -0 m ) = l-u*(l-s)-?;*(l-s 2 ) if \8 S - 0*\D O l < R* (10) 
= otherwise, 
where s = y/l - {\0 S - 0* \D l/R*) 2 , 

and where 0* is the angular position of the star's center (recall that we have chosen the origin to be at 
the center of the lens, the planet). The limb darkening parameters it* and u* are generally wavelength 
dependent. Typical values are u* = 0.8, v* = —0.225 in red bands, it* = 0.99, = —0.17 in blue bands, 
and m* = 0.93, = —0.23 in intermediate yellow bands (Cox 2000). Clearly, as far as the gross features of 
the transit are concerned, it is the ratio R*/Rq that is important, not the absolute size of i?* (other than a 
rescaling in the overall duration of the transit). 

In summary, we have: three lensing parameters, e, B, and Rq/H; two transit parameters, #tr.imp., and 
7tr.imp.; and three stellar parameters, R*/Rq, u*, and v*. In the case of a spherical planet the total number 
of parameters is reduced by two because the lensing paramter e and the transit parameter 7t r .imp. are n °t 
needed. We will not perform an exhaustive study of the full parameter space in this paper, but will be 
content with a mostly qualitative (with a few exceptions) description of the dependence of the transit light 
curve on these parameters. Some of these parameters arc likely degenerate. For instance, changing $tr.imp. 



3 The parameter fltr.imp. i s related to the planet orbital inclination i by 6*tr.imp. = [Dls/D ol]cosi. 
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means the transit is sampling a different part of the stellar profile, which might be mimicked by a different 
limb darkening law (e.g. a different and v*). Finally, note that for extrasolar systems, -Dol 3> -Dls and 
so -Dol drops out of the lensing equation (e.g. equation (||)), and appears only as an unobservable overall 
scaling of angular separations. 



4. Strong Lensing and Caustic Structure 

4.1. Condition for Strong Lensing 

The first question we would like to address is when strong lensing — i.e. the existence of caustics or 
multiple images (which may be unresolved) — occurs. The existence of caustics could lead to significant 
modifications in transit light curves. 

Caustics in the source plane, or critical curves in the image plane, can be obtained by solving for 
divergent magnification, A^ 1 = (equation (0)). In the spherical case, with e — 0, this is straightforward. 
Recall that images must have u > Ro/Dol to be visible, otherwise they are blocked (occulted or absorbed) 
by the planet (equation (^J)). Therefore, the condition for strong atmospheric lensing by a spherical planet 
is 

1 + 2V>(u) + 4>{uf + u 2 tP{u)^(u) + u 2 ip{u) = [1 + Tp{u)][l + i)(u) + u 2 ip{u)} < 0, u = Rq/D l- (11) 

This guarantees that some images with u > Rq/Dql will have A^ 1 — 0. Using the results in §||, the above 
condition imposes a relation between only two parameters: B and Rq/H. 

Since e is small, the condition for strong lensing by an elliptic atmosphere will not be too different from 
the spherical case described above. We estimate it using results proved in Appendix B (equation (|3^) and 
the following paragraph), replacing the above condition with the following: 

1 + 2i/>(u) + (1 - e 2 )t/j(uj 2 + (1 - e a )uV(u)V!(u) + (1 - e)u 2 ^(u) < 0, u = R /D OL . (12) 

Using equation this imposes a condition on three parameters, e, B and Rq/H, 

l -fi{i B - f<0 - < 13 > 

where we have used H/Ro <C 1 and eCl. This is a main result of our paper. 

The above condition for strong lensing is depicted in Figure |3[ We show in the same figure the relevant 
parameters for several known planets. We assume a = 1.2cm 3 /g for all of them. Mars (M) has i?o ~ 3400 
km, D oh D hs /{DoL + D LS ) - 1-5 AU, p Q - 1.5 x 10- 6 g/cm 3 and H - 8.7 km (Jones 1999 p. 314). Jupiter 
(J) has R ~ 71000 km, D OL D LS /(D OL + D LS ) - 5.2 AU, p ~ 3.5 x 10~ 5 g/cm 3 and H - 21.7 km (Jones 
1999 p. 341). The extrasolar planet HD 209458b (HD) has i? - 10 5 km, D lD ls /{D 0L + D LS ) - 0.05 AU, 
pa - 1 - 6 x 10~ 6 g/cm 3 and H ~ 500 - 700 km (Charbonneau et al. 2000, Henry et al. 2000, Mazeh et al. 
2000, Jha et al. 2000, Burrows et al. 2000). The density po is obtained by computing the density at which 
the optical depth for Rayleigh scattering is unity (at A = 6500 A; see §||), or in the case of Mars, po is the 
density at the base of its atmosphere, just above its solid surface. For a gaseous planet or a rocky planet with 
a thick atmosphere, po (and Rq) is wavelength-dependent. For example, our own solar system planets have 
much stronger molecular absorption in the UV and IR compared to the optical. Strong absorption bands 
make po smaller (where, again, po is the density and Rq the radius at which the optical depth = 1), in some 
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cases to the point where atmospheric lensing is no longer strong (recall that B oc po). This is illustrated by 
the point J' where we show Jupiter observed at a waveband where the absorption cross-section is 600 times 
larger e.g. Rayleigh scattering observed at A = 1300 A instead of 6500 A. The main effect is a dramatic 
decrease in B. The increase in Rq/H is very small, because the radius varies only logarithmically with the 
density. We will discuss the color dependence of lensing in more detail in §^|. 

There is clearly a wide range of possibilities in the two-dimensional parameter space of B and Ro/H. J 
and HD, which have a similar mass, occupy very different parts of the diagram primarily because they are 
situated at very different distances from their parent stars. Higher temperature in the latter leads to a much 
larger atmospheric scale height. To investigate how variables such as stellar luminosity, albedo, planet mass 
and composition control the lensing behavior is outside the scope of this paper. In fact, the radial velocity 
detections of extrasolar planets have taught us that planet orbital characteristics can be quite different from 
expectations based on the solar system. We therefore adopt the philosophy that the parameters spanned by 
the above three examples are all possible and interesting, and our goal is to understand the importance of 
atmospheric lensing under these conditions. We will briefly discuss the physical motivations for the choice 
of some parameters in §M. 



4.2. Caustic Structure and Magnification 

The lower panel of Figure |i] shows the caustic (solid line) and the critical curve (dotted line) for a 
fiducial model with projected oblateness e = 0.05 denoted by the solid square in Figure^ (see §[| for details). 
A point source situated within the caustic produces 4 images 

The upper panel of Figure ^] shows the magnification as a function of source position for a point source 
situated on the x-axis. This figure shows that lensing has two effects. First, it suppresses the flux originating 
from source positions just outside Rq (uDol > Rq)- Second, it brings in additional photons from source 
positions behind i? {uDol < Ro) — photons which would otherwise be blocked (occulted or absorbed) by the 
planet. Which effect dominates can be calculated by integrating over the star (regarded here as a collection 
of point sources) and depends on details of the magnification profile and the limb-darkening profile (i.e. 
the lensing parameters, transit parameters, and stellar parameters described in fj3|). We can derive a simple 
result from equation (Q) for the limiting case of a completely flat and constant /. Taking / out of the integral 
in equation (^), F can be rewritten as / J d 2 9iW(0i), which is identical to the result if no lensing takes 
place. It is not hard to see that the same conclusion follows if I is constant within some region, and if the 
planet is well within this region. In other words, if there is strictly no limb darkening, there will be no net 
gain or loss of stellar flux from atmospheric lensing, provided the projected planet is well within the stellar 
disk. 

Realistic stellar profiles are never exactly constant. Refraction therefore generally modifies the dimming 
or deficit of the stellar flux during a transit, especially during ingress and egress. The sign of the modification 
depends on the exact transit and lensing parameters. The size of the modification can be crudely estimated 
from the ratio of the projected area of the atmosphere to that of the planet: ~J7/ Rq, if strong lensing occurs. 



4 It is interesting to note that lensing by an elliptic potential or density profile generally produces two sets of caustic curves 
(e.g. Schneider, Ehlers & Falco 1992), whereas we see only one here. The origin can be traced to the fact that tp(u) in equation 
(H), which is proportional to the deflection angle, is monotonic. We do not expect small local deviations from the exponential 
density profile to change this conclusion, because ip(u) involves an integral over many layers of the atmosphere. 
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4.3. Application to Solar System Occultations 

Figure ^ also represents an occultation of a point-source background star by a solar system planet. The 
dotted line shows what one would expect for an atmosphereless planet (ignoring diffraction): the stellar 
intensity is constant and vanishes instantaneously when the star passes behind the planet as viewed by the 
observer. The dashed line shows the magnification from equation (Q) when refraction by the atmosphere 
is significant. The magnification curve can be understood qualitatively as follows. (1) The diminution of 
starlight during ingress and egress is due to atmospheric differential refraction (not by absorption which 
occurs at much lower depths in this model) , which causes stellar light rays to diverge away from the planet- 
observer line-of-sight. This is a consequence of the radial density gradient in the atmosphere. (2) When 
the star is behind the opaque part of the planet, the occultation light curve is not zero, because some of 
the stellar light rays are refracted into the observer's line-of-sight. (3) When the point source is near the 
geometric center of the planet, the stellar flux is symmetrically refractively focused causing an increase in 
brightness. The solid line in Figure || shows the occultation light curve that includes the flux from the planet 
from reflected sunlight if the occulted star has an intensity 20% of the planet. The solid curve shows that the 
planet's reflected light dominates the light curve during the stellar occultation. Hence the minimum observed 
brightness corresponds to the planet alone and the full extent of refractive focusing and defocusing is not 
observed for bright solar system planet occultations of background stars. The magnification peaks have been 
observed for most solar system planets (and with the multiple images even resolved in the case of Saturn; 
see Nicholson, McGhee, & French, f995). Note that the magnification peaks are generally reduced due to 
atmospheric absorption or scattering and to the finite size of the star. Note also that if the star does not 
move along the axes of symmetry of the solar system planet, the two magnification peaks will generally be of 
different heights and widths. Also, the two peaks will merge into one if the planet's oblateness is sufficiently 
small. Refractive occultations of all of the solar system planets (except Mercury) have been observed and 
in some cases well studied (see Elliot & Olkin 1996 and references therein). Models of them can provide 
temperature (T), pressure (P), and density (p) as a function of atmospheric depth. 



5. The Transit Light Curve 

5.1. Exploring the R Q /H - B Plane 

Figure |^ shows the transit light curve with or without lensing for a fiducial model denoted by the solid 
square in Figure |: B = 40.3, R /H = 117.3, e = 0.05, together with P /P» = 0.084, R*/y/GM Q /D hS = 
233.66 minutes, it* = 0.8 and v* = —0.225, (9 tr .i mp . = 2Pq/Pol and 7tr.imp. = 45°. The last six parameters 
define the overall light deficit and rough duration of transit, the limb darkening profile, and the geometry of 
the transit. We emphasize that a different value for R*/y/ GMq / Pls (§[§) can be easily accommodated by 
rescaling the time axis of all light curves shown below. 

To better isolate the effect of lensing, we define and show in Figure || the following quantity (solid line): 

f(t) = [AF lcns (t) - AP nolcns (i)]/AP/° a s x -, (14) 

where AF(t) = 1 — F(t) represents the light deficit, and the subscripts denote whether or not lensing is 
taken into account (with all other parameters fixed). We choose to normalize by the maximum deficit, which 
can be estimated roughly by (Pq/P*) 2 -0 The quantity / is therefore a fractional difference. This gives us 



5 We could normalize by the maximum deficit in the nonlensed case, but it would only modify f(t) to higher order (0(/(t) 2 )). 
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a nice way to estimate the absolute size of the deficit for a transit with any given Rq/R*: it is roughly 
(i?o/-R*) 2 (l + /)• A s we have discussed in §[|, / can be estimated by ~ H/Rq, which is around 1% for 
our fiducial model. This is in rough agreement with the results of our numerical integration. The absolute 
change in flux due to lcnsing is therefore 1% x (R /R r ) 2 , which for the example shown with R /R* ~ 0.1, 
is 10 -4 — a potentially significant effect. 

The solid line in Figure ^ shows the same quantity / for a transiting planet with the same H/Rq, 
but a smaller B (= 14.9) — denoted by an open square in Figure ||. The overall magnitude of / is a little 
smaller but not too dissimilar from the fiducial model above. The sign of the lensing effect as the transit 
progresses, however, is quite different. As we have previously discussed, two opposing lensing effects are 
competing. One is the suppression of flux from source positions just above Rq. The other is the addition of 
photons from source positions in the planet's shadow. For the model with a larger B, the first dominates 
over the second effect at all times. For the model here with a smaller B, the second effect dominates, at least 
temporarily towards the end of ingress or beginning of egress. At mid-transit, however, the first effect still 
dominates. Interestingly, a search through a whole range of limb-darkening parameters reveals that lensing 
almost always causes a net suppression of flux at mid-transit, except when the stellar profile is very spikey 
at the center. 

When B is decreased further so that caustics no longer exist, the first effect, suppression of flux, will 
dominate. This is illustrated by the dotted line in Figure fj], which is the model denoted by a cross in Figure 
|3| and is chosen to resemble HD 209458b (see § [4.1[ ) f\ Clearly, atmospheric lensing is weak for the close-in 



extrasolar giant planets like HD 209458b. This is no surprise since, as we have pointed out in §4.1 and Figure 
||, there are no caustics in the case of HD 209458b. This is due to the low B which is mainly caused by a 
small lens-source distance because B is proportional to Dls ■ A second, less direct effect of the small Dls is 
that the planet's temperature is hotter and thus H is higher, further decreasing B. However, it is interesting 
that even though the combination of parameters does not allow the existence of caustics in the case of HD 
209458b, refraction nonetheless could modify the light curve at a (fractional) level of ~ 5 x 10" 4 (depending 
on the actual atmospheric structure; see §[?] for details). Note also that because Dls = 0.05 AU for this 
close-in gas giant (compared to Dls = 1 AU we adopt in previous cases); the shorter orbital radius leads 
to a higher velocity (the parameter w; see §[|), hence a shorter transit. Finally, the dashed line in Figure R 
f(t) shows the model shown by a triangle in Figure ^|. The behavior is similar to the solid line of Figure ^ 
except that the overall effect is weaker. This is due to a larger Rq/H, ~ 700, from which one expects / of 
the order of ~ 1.5 x 10~ 3 , in agreement with what we find. 



5.2. The Lensing Signal 

The previous section shows that lensing can modify the transit light curve to an extent that is potentially 
detectable in some cases. If so, high precision photometric observations of extrasolar planet transits can 
provide planet atmosphere parameters such as scale height and density, both of which affect the lensing 
behavior. What is not clear, however, is whether a transit light curve from a planet with lensing can be 
distinguished from a transit light curve from a planet with slightly different parameters, but with no lensing. 
This is what we investigate here. 



6 The projected ellipticity is chosen to be 10 — 3 , although its precise magnitude does not affect our conclusions here. The 
small oblateness is expected due to the tidal locking of rotational and orbital motions, given the proximity of HD 209458b to 
its companion star (see Hui & Seager, in preparation). 
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As shown in the last section, the modification introduced by atmospheric lensing is qualitatively quite 
different for models with or without caustics. For models with caustics, such as those shown in Figure ^, 
f(t) experiences a significant dip in the light deficit around ingress and egress (the same holds for solid and 
dashed lines in Figure 0). This dip is caused by the additional photons brought in from the planet's shadow 
due to lensing which causes the light deficit to drop temporarily. In other words, the inexorable dimming of 
star light at ingress is temporarily slowed down by the additional photons that are refracted into view. This 
effect is much weakened if caustics do not exist, as in the case of the dotted line in Figure ^, which is our 
model for HD 209458b. This dip, or slight reduction in the rate of dimming, introduces a distinct shape to 
the transit light curve during ingress and egress, which one might hope to observe. 

The dotted line in Figure ^ shows the following fractional difference: 

f(t) = [AF lcns (t, P) - AF nolcns (t, P')]/AF^(P), (15) 

where the extra argument P denotes the whole set of parameters that determine the transit characteristics. 
We use P' for the second term (which does not take into account lensing) to emphasize that the second term 
has a different set of parameters from the first term. The question is how small one can make / by choosing 
P' appropriately. If f(t) is too small to be detectable, then for practical purposes, one cannot tell from the 
observed lightcurve whether or not refraction plays a role, and therefore cannot determine useful parameters 
from refraction effects. 

The dotted line in Figure || shows the result of our effort of searching for the appropriate AF no i ons (t, P') 
that would minimize / (t) for our fiducial model (solid square in Figure |J) . We systematically vary a host of 
parameters, including Rq, u», v„, w, #tr.imp. an d Ttr.imp.- The minimum fractional difference we can come 
up with is a model that brings f(t) down to about 10~ 3 during mid-transit, but still almost 5 x 10~ 3 at 
ingress and egress. It should be emphasized, however, that it is conceivable a different limb darkening law 



from the one we have adopted (equation (10)) could be chosen to further decrease /. Nonetheless, within 
the range of simple models we consider, it appears that the lensing signal (i.e. the caustic signal) cannot be 
completely masked by a clever choice of parameters. A discrepancy in / of 5 x 10 -3 is potentially observable; 
this fractional difference of 5 x 10~ 3 means, for a transit where Rq/R% ~ 10~ 2 say, the absolute difference 
in the light deficit would be ~ 0.5 x 10~ 4 . We will discuss other nonlcnsing effects that might be of similar 
magnitude in 

We find that the same conclusion holds for other cases with caustics, such as those shown as solid and 
dashed lines in Figure |7[ For cases without caustics, on the other hand, such as that shown with a dotted line 
in Figure J?], we find that a suitable choice of Rq for the second term in equation ( |l5| ) is generally sufficient 
to make / quite small and undetectable. 

The time derivative dF/dt provides another way to illustrate the effect of refraction. Figure || shows 
the time derivative dF/dt for our fiducial model, with (solid line) or without (dotted line) lensing. The 
bottom panel focuses on the part of ingress where dF/dt reaches its minimum value before increasing. This 
corresponds roughly to when half of the planet has crossed the stellar limb. What is interesting is that the 
presence of caustics introduces a distinct enhancement in dF/dt (solid over dotted curve) around t = —230 
min. This is due to the additional photons brought about by the caustic. From Figure ||, we can see that 
the caustic has a size of about 0.2i?o/-DoL- Dividing it by w, the angular velocity, we obtain ~ 5 minutes, in 
rough agreement with the duration of the enhancement seen in Figure |^. Since the size of the caustic scales 
with e, this offers a direct way of measuring the projected oblateness of the planet. However, it appears 
that one needs to be able to detect differences in dF/dt at the level of at least one part per million per 
minute, which is likely difficult to achieve in the near future. Therefore, while the derivative dF/dt offers 
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an interesting way to look at the effect of lensing, the integral fractional differences / or / might be a more 
practical probe. 



5.3. The Oblateness Signal 

The projected oblateness or ellipticity of a planet can affect the transit light curve. Here we explore the 
fractional difference 

f e (t) = [AF(t, e - O.f ) - AF(t, e = 0)]/AF max (e = O.f ) (16) 

to determine the magnitude of the effect. Here AF can either have lensing taken into account or not because 
we would like to investigate how oblateness affects the light curve irrespective of whether or not refraction 
is present. Figure ^ shows f e for the fiducial model (solid square in Figure ||[). The solid line denotes a case 
with lensing, and the dotted line without. It is interesting that oblateness introduces a modification to the 
light curve that is quite similar in the two cases. 

Next we pursue the same exercise as before and ask if parameters P' can be found such that the following 
can be minimized: 

f e (t) = [AF(t, e = 0.1, P) - AF(t, e = 0, P')]/AF™( £ = 0.1, P), (17) 

where P denotes the parameters for our fiducial model. Figure [l^ shows the results. We cannot find 
combinations of parameters P' to reduce f e below I0~ 3 at ingress and egress, suggesting that we might be 
able to distinguish light curves caused by ellipsoidal versus spherical planets. The P' that minimizes f e , as 
shown in Figure [l^, turns out to correspond to a spherical planet with slightly larger Rq, and essentially the 
same area as the ellipse in the fiducial model. In other words, the curves in Figure |l0| basically show the 
fractional difference between the light deficit of a spherical planet and an oblate planet (e = 0.1) with the 
same area (with or without lensing). This is a detectable signature. 

How large is the lensing signature compared to the oblateness signature for the same stellar parameters? 
For our fiducial planet model, a comparison between Fig. ^| and (dotted line of) Fig. [T(] shows that the 
two effects are of a comparable magnitude. Note, however, the oblateness signature can persist even if the 
parameters are not right for strong lensing. 

The reader might have noticed that some of the previous light curves are not symmetric around t = 
(which is chosen to be the point of mid-transit or minimum flux). This is due to the fact that the major/minor 
axis of the planet is misaligned with the direction of motion (see Figure ||) i.e. only 7t r .imp. = or 7t r .imp. = 
90° would produce a light curve that is symmetric around t = 0. In the presence of a general misalignment, 
the asymmetry provides another useful diagnostic of oblateness because a spherical planet can only produce 



symmetric lightcurves. We show in Figure 1 1 the following fractional difference, 



/ asym (t) = [AF(t) - AF(-t)]/ AF max . (18) 

The solid line denotes the case for e = 0.1, and the dashed line for e = 0.05, both for 7tr.imp. = 45°. The 
dotted line shows the same for 7t r .imp. = 0, verifying that the asymmetry vanishes for exact alignment of 
major (or minor) axis with the direction of motion. The lower panel shows f aS ym(t) where all quantities are 
evaluated with lensing, while the upper panel shows the same without lensing. The sets of curves are almost 
identical. This means the degree of asymmetry in the light curve is determined by the size of the oblateness 
and angle 7tr.imp. alone, and is insensitive to lensing. 
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6. Absorption and Color Dependence 

6.1. More Accurate Modeling of Absorption and its Effects on Lensing 

So far we have been using a crude step function model of absorption that turns abruptly on and off at 
u = Rq/Dol- Here we use absorption to refer to extinction by either absorption or scattering. Taking into 
account absorption in a more realistic way is in principle straightforward. For each image position, there is 
an associated optical depth: 

/oo 
a— P —dz (19) 
-oo P m H 



Po 



■a^/2TruD OL H exp[-( Uj D OL - Rq)/H], 



where u is defined in terms of Bi as in equation (|^), and we have used the condition that uDol/H 1. 
The symbol fimn denotes the mean molecular weight. To compute the light curve with absorption, W(9i) 
should be replaced by e~ T ^ I ' ) in equation (0). 

The absorption cross-section a can vary significantly depending on the composition of the planet at- 
mosphere and the wavelength of interest, and is essentially infinite at the surface of a rocky planet. To 
demonstrate concretely the effect of more accurate, gradual absorption, we adopt the cross-section appro- 
priate for Rayleigh scattering (Jackson 1999): 

a = l(T 27 cm 2 (5000A/A) 4 , (20) 

where A is the wavelength of interest. 

We recompute predictions for the fiducial model (solid square in Figure |^) using the model outlined 
above. In Figure |l2|, we show the following fractional difference, 

f a (t) = [AF ah (t) - AF(t)]/A^ ax (21) 

where AF a b(£) represents the deficit computed when a gradual change in absorption is used, while AF(t) 
denotes the same computed with a step function approximation for W (equation (^|)). The upper solid and 
dotted lines in Figure [l^ give the above quantity with or without lensing. Gradual absorption produces 
light curves that are different from a step function absorption/occultation at about the 1% level. The above 
uses the same planet radius i?o for AF a b and AF. This difference can be much diminished, however, if one 
chooses an appropriate i?o for the step function model. The lower solid and dotted lines in Figure |l2| are 
the minimized fractional difference: 

f a (t) = [AF ab (i, R ) - AF(t, R' )]/AF^(R Q ), (22) 

for cases with and without lensing respectively, where we tune R' to make f a small. The small value of f a (t) 



in Figure 12 shows that our step function absorption model is actually a good approximation to reality — a 
suitable Rq can always be found such that it approximates a gradual absorption model to high accuracy. This 
works primarily because e~ T does behave almost like a step function, due to the fact that p and therefore r 
varies exponentially with radius. 

With the above results we can proceed to discuss the color dependence of a transit light curve. There are 
three main factors. First, the stellar profile changes with wavelength, generally flatter as one considers redder 
wavebands (see §0). Second, the refractive coefficient a varies with wavelength. The variation depends on 
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atmosphere composition, and some examples are given in Table |l| Third, the absorption cross-section a also 
varies with wavelength. In fact, for Rayleigh scattering for instance (equation (|20|) ), absorption varies much 
more strongly with wavelength compared to both the index of refraction and limb darkening. This can have a 
dramatic consequence for the existence of caustics. Consider a model like the open square in the B — Rq/H 
plane as shown in Figure [| Recall that B = [2ap /H] [-Dql-DlsA-Dql + -Dls)] (equation (^)). As one 



considers shorter wavelengths or bluer colors, a becomes larger, increasing the optical depth (equation (19)). 
In our step function model, this is equivalent to lowering po (defined to be the density at which r = 1), 
in other words lowering B^ For sufficiently blue colors, the model would shift down from the open square 
in Figure || and cross the strong lensing threshold (the solid line), erasing the caustics, and making the 



lensing signal much weaker (as discussed in §5.2). The same also holds true at wavebands where other 
kinds of absorption are important, such as molecular electronic absorption bands in the UV, and rotational- 
vibrational bands in the IR of species such as EbO, CO2, and CH4. This is a well known result for solar 
system planet occultations of distant stars where lensing effects are strong in the optical but nonexistent in 
the UV and within strong absorption bands in the IR. 

In summary, for extrasolar planet transits where caustics exist, we expect the lensing signal to disappear, 
or at least weaken, at wavebands with high extinction. 



6.2. Planetary Atmosphere Transmission Spectrum and Stellar Spectrum 

Extrasolar planet transit transmission spectra have been described in several papers for both the close-in 
extrasolar giant planets (Seager & Sasselov 2000; Brown, 2001; Hubbard et al. 2001) and Earth-like planets 
(Schneider 1994; Webb & Wormleaton 2001). For parameters with strong lensing (Figure 2) the transmission 
spectrum could be significantly affected by atmospheric refraction. At wavebands corresponding to the 



transparent continuum, strong lensing can lower the flux (recall from §5.1 that lensing generally causes a 
net suppression of flux during most of the transit), while as we have discussed in the last subsection, at 
wavebands corresponding to strong planetary absorption bands, lensing has negligible effects. The contrast 
between continuum and absorption lines is therefore decreased, and effective line-strengths are therefore 
altered. 

The magnitude of this effect can be estimated as follows. Atmospheric lensing effect generally introduces 
a decrease in absolute flux that is approximately (Rq/R*) 2 x (H/Ro) for planets with caustics. Strong 
planetary atmosphere absorption lines on the other hand cause a change in absolute flux that is approximately 
(i?o/-R*) 2 x x x (H/Rq), where x ~ a few. Therefore, in a planetary transmission spectrum, the effective 
line strength is reduced from x x (H/Rq) to ~ (x — 1) X (H/Rq). The exact size of this effect depends on 
details of the atmospheric structure, composition as well as the limb-darkening function. 

As noted in our Figure ||, and commented on previously (Seager & Sasselov 2000; Hubbard et al. 2001; 
Brown 2001) the transmission spectra of the close-in extra-solar giant planets (EGPs) — including the only 
known transiting EGP HD 209458b — will be little affected by refraction because caustics do not exist due 



to the small planet-star distance (see §5.1) 



7 Changing a also changes our definition of Rq which is where r ~ 1. While this does affect the transit light curve somewhat, 
its effect on lensing is smaller compared with that due to po- This is because varying <r typically changes Rq by a few scale 
heights H, and that represents a small change to Rq/H (since Rq 3> H), which is the other parameter that controls the lensing 
behavior. 
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Finally, we end this section by briefly commenting on the stellar spectrum. During a planetary transit, 
the stellar spectrum changes with time, a phenomenon referred to as a spectroscopic transit. There are at 
least two different effects. The first has been discussed and observed by Queloz et al. (2000). The planet 
blocks different parts of the rotating stellar disk as the transit progresses, causing red-shifting or blue-shifting 
of the stellar lines during the transit. The second effect arises from the fact that different parts of spectral 
lines and different lines are formed at different depths in the stellar atmosphere. As the planet transits 
different parts of the limb-darkened stellar disk, the stellar line shapes and strengths changes. Atmospheric 
lensing introduces modifications to both of the above effects, because flux from different parts of the star is 
magnified compared to the nonlensed case. 



7. Discussion 

Our findings are summarized as follows. 

• The importance of atmospheric lensing is mainly controlled by two parameters: Rq/H (the ratio of 
planet radius to atmospheric scale height) and B = 2apoD^s/ H (product of the refractive coefficient, 
atmospheric density and star-planet separation divided by the scale height) .^] The condition for strong 
lensing — the existence of caustics — is described by equation ( |l3| ) and depicted in Figure ||[ 

• Strong lensing generally introduces a fractional change in the light deficit of the order of H/R during a 
planetary transit. We choose to discuss changes induced by lensing in terms of fractional changes in the 
light deficit (f(t) as defined in equation ( fhij ) and shown in e.g. Figure^) so that the absolute change 
in observed flux can be easily estimated for any planet-to-star size ratio.^] It works as follows. The 
drop in flux during a nonlensed transit is about (Rq/R*) 2 where i?* is the radius of the star. Strong 
lensing introduces an additional absolute change in flux that is therefore about (Ro/R*) 2 x (H/Rq). 

• Lensing produces a characteristic slowing down in the dimming of the star light at ingress (reversed 
at egress). This is due to the additional photons refracted into view from the planet's shadow. The 
light curve of a lensed transit can be mimicked to some extent by one of a nonlensed transit, if 
parameters for the latter are appropriately tuned. We find, however, that the difference between the 
two can remain significant especially during ingress and egress. For example, the dotted line in Figure 
^ shows a fractional difference in the light deficit of ~ 0.5(H/Rq) at ingress or egress. When this 
difference is larger than the observational uncertainties, we should be able to constrain the lensing 
signal parameters such as the atmospheric scale height, density and refractive coefficient (parameters 
such as i?o and -Dls can be learned from other features of the light curve or other observations). The 
scale height and refractive coefficient will in turn give us useful information on the temperature and 
chemical composition of the planet. 

• The strength of lensing is expected to vary significantly with color. The primary reason is the variation 
of absorption and scattering with wavelength. In wavebands where absorption or scattering is signif- 
icant, caustics are effectively hidden by extinction. Therefore, the optimal wavebands for detecting 
atmospheric lensing is towards the red or near infrared, where the inevitable Rayleigh scattering is less 



In the solar system case, Dls is replaced by -Dql, see 




The observed stellar flux here is always normalized by its pre- or post-transit value. 
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strong, but away from strong absorption bands such as gaseous H2O or CH4. Optically thick clouds 
(such as enstatite in hotter planets or H2O or CH4 ice in cooler planets) might contribute significant 
opacity. See Hubbard et al. (2001) for a detailed treatment of HD 209458b. 

• Several conditions are therefore advantageous for detecting the lensing signal: 1) observations at wave- 
bands where absorption and scattering is weak; 2) a high temperature of the planet atmosphere due 
to a nearby hot star, which raises the atmospheric scale height (thereby increasing H/Rq); 3) Dls 
large enough to keep B = 2apoDi J s/H above the strong-lensing threshold. Our fiducial model (solid 
square in Figure |J) provides an example where lensing induces a fractional change in the light deficit 
of around 1%. For a planet that has an area about 1% of the star, this implies lensing causes an 
absolute change of about 10 -4 in the normalized flux — an observationally relevant effect. This model 
can be realized by a gaseous giant (similar to Jupiter in mass, but with a slightly inflated radius due 
to high temperature) about an AU away from an A star, with observations done in red wavebands 
(A ~ 10 4 A), but away from possible strong absorption bands such as water. We should emphasize 
there is considerable uncertainty in the size of such a planet, because the radius of a hot gaseous giant 
varies significantly with time (Burrows et al. 2000; see also Murray et al. 1998 and Lin et al. 1998 on 
change in orbital characteristics over time). 

• Oblateness of the planet induces an asymmetry to the transit light curve (about the point of minimum 
flux), which vanishes only when the semi-major or semi-minor axis of the planet is exactly aligned with 
the direction of relative motion. The asymmetry is about 0.5%, as measured in fractional light deficit 
(equation ( |l8| ) and Figure (|TT|) ) , for a projected oblateness or ellipticity of e = 0.1, and an angle of 45° 
between the semi-major axis and direction of motion. The absolute difference in flux between ingress 
and egress is therefore about 0.5% x (Rq/R*) 2 , and is for example 5 x 10" 4 for a planet that is a tenth 
the size of the star. 

• For HD 209458b, the only currently known transiting extrasolar planet, caustics and therefore strong 
lensing is absent, because of the very small lens-source separation (and the resulting large scale height 
due to a high temperature). Its oblateness is expected to be small, < 10~ 3 , due to the tidal locking of 
its rotation to orbital motion. 

In this paper, we have focused on the effects of atmospheric lensing and planet oblateness on the transit 
light curve. It is worthwhile to briefly list other possible "secondary" fluctuations to the light curve — 
variations in the observed flux other than those due to the standard spherical occultation of a star with a 
smooth limb darkening profile. They can be divided into three categories. 

The first category is due to close companions of the planet. Rings and moons, if present, induce 
an asymmetry to the light curve that may be confused with oblateness. They also introduce additional 
parameters that control the light curve, making it perhaps more difficult to discern the lensing signal. 
Taking Jupiter as an example, the largest moon Ganymede has a radius of 2635 km, about 3.7% the size of 
Jupiter. The same moon of a Jupiter-like planet orbiting another star would induce a fractional change of 
about 10~ 3 in the light deficit. This is about 5 — 10 times smaller than the lensing or oblateness effects we 
found for our most optimistic models (Figures |^ and [Tl| ). If the moon is sufficiently far from the planet, it 
might also introduce distinct signatures at ingress or egress (e.g. Sartoretti & Schneider, 1998) that can be 
disentangled from oblateness or lensing effects. An opaque ring 1/2 the size of Saturn's rings for a gaseous 
giant with R p = 1.4 Rj orbiting a sun- like star would deepen the planet transit light curve by 0.05% - 1.2% 
depending on the inclination of the ring (Seager & Hui 2000). Ring effects might be isolated by looking for 
distinct shapes in the transit light curve (see Brown et al. 2001), and possibly rescattered star light. 



-17- 



The second category is due to the planet itself. The planet atmosphere can have persistent or transient 
disturbances that modify the transit light curve. Such disturbances have to be large in extent to change the 
light curve significantly. For instance, small uniformly distributed clouds would not introduce an asymmetry 
to the light curve. Something like the Great Red Spot on Jupiter could conceivably be large enough. But 
such disturbances cannot introduce a fractional change in the light deficit that is larger than the ratio of 
the projected atmospheric area to total planetary area i.e. ~ H/Rq. It is important to emphasize that 
perturbations to the atmosphere that only change the effective radius of the planet Rq are not sufficient to 



wash out the lcnsing signal; as we have shown in i 5.2 , just tuning Rq is not enough to confuse a lensed transit 
with a nonlensed transit that has different parameters. Hubbard et al. (2001) pointed out that Rayleigh 
scattering, in addition to causing extinction, also produces a glow around the planet, which changes its 
effective size (fractional variation in size about 1% between different wavelengths). It would be interesting 
to explore how both Rayleigh and condensate rescattered stellar photons affect the detailed shape of the 
transit light curve (i.e. aside from a simple change in i?o)- Lastly, we have also ignored diffraction in this 
paper. Diffraction is likely unimportant for most extrasolar cases of interest, for two reasons. First, the 
relevant parameter combination (R/H^^BDol/Dls is greater than unity (see Fig. ||), and is therefore 
not favorable for diffraction (see e.g. French & Gierasch 1976, Elliot et al. 1975). Second, the observed flux 
comes from a sum over incoherent sources distributed over the stellar surface. 

The third category of secondary fluctuations is due to complications in the star. Stellar oscillations 
are expected to produce absolute changes in the flux at a 10 -5 level on the time-scale of a transit (see e.g. 



website for the proposed Kepler mission ittp : //www . kepler . arc . nasa . gov/ ) . Also, realistic stellar profiles 



might have bumps and wiggles on top of the smooth profile we have assumed (§10j). For instance, the light 
curve would be modified if a planet happens to transit over a star spot (e.g. Seager & Hui 2000). The largest 
effect is obtained if the spot has a size that is comparable to the planet, and the planet happens to overlap 
completely with the spot during the transit. At the point where the planet and spot coincide, the flux would 
basically return to unity, its pre-transit value. For example, an Earth-sized planet which happens to cross 
over an Earth-sized star spot would produce this behavior, as illustrated in the upper panel of Figure [l^. 
This kind of situation is probably rare, however, particularly if the planet considered is a gaseous giant i.e. 
much larger than typical star spot (lower panel of Figure |l3| ). Then, the absolute change in flux in such 
coincident transits will be roughly the ratio of the spot area to star area, < 10 -5 borrowing from the example 
of the Sun. Note also that while lensing or oblateness effects are most easily recognized at ingress or egress, 
changes to the light curve due to imperfections in the stellar profile can occur throughout the transit. More 
detailed studies would be required to see if spots produce different signatures from lensing or oblateness, if 
the planet crosses them at ingress or egress. Note also that effects due to star spots are likely not repeatable 
over many planetary periods. 

It is also worth emphasizing that even in cases where the lensing or oblateness signal is weak and 
therefore not easily isolated, it could still act as an important contaminant in confusing other signals one 
might be interested in, such as detection of moons and planetary rings. From the above discussions, it 
is clear future high precision measurements of extrasolar planetary transits will present a very interesting 
challenge — there are several sources of small secondary fluctuations (change in flux of 10~ 4 or less) which are 
within observational reach and which would require some effort to disentangle. The rewards of such efforts 
will be substantial — from detection of moons and rings, to measuring the (projected) planet oblateness and 
therefore constraining its rotational period; from studies of the stellar atmosphere to studies of the planetary 
atmosphere, constraining for instance its scale height, temperature, density and composition. 
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Appendix A 

Here, we would like to motivate the density profile we adopt in this paper. We begin by the following 
statement of hydrostatic equilibrium of the atmosphere: 

VP = -pV^ cff . (23) 

where P is the pressure, and <p e g. is an effective potential which is the sum of the gravitational potential (j) 
and a potential to take into account rotation: 

0cff. =0-^ 2 (x' 2 +2 / 2 ) (24) 

where z' is the axis of rotation, x' and y' are the perpendicular axes (we will use x[, i = 3, 1,2 to denote 
them) and lo is the angular speed. For an isothermal atmosphere where P oc p, the above implies 

p oc exp[-(/> cff .] (25) 

We note that planetary atmospheres are only approximately isothermal. Deviations such as temperature 
inversion are known to exist (e.g. Jones 1999). 

The question then reduces to what kind of gravitational potential one expects for a general rotating 
figure of equilibrium. For bodies that are not too aspherical, we can expand the gravitational potential using 
the Legendre polynomials (see e.g. Danby 1962, Chandrasekhar 1969) 



l--^-(3cos 2 tf-l)-...) (26) 



where M p is the planet mass, J 2 is some constant coefficient which is presumably small, 9 is the angle between 
the radial vector and the z' axis, and r' = \J x' 2 + y' 2 + z' 2 . It is sufficient to illustrate our argument using 
only the first term - it is possible to generalize to include higher order terms in the potential expansion. 

We are interested in the form of (j) c s. when r' = R(l + i]) with i)<l, where R is approximately the 
planet radius. We are also interested in cases where the angular speed u> is in some sense small — the 
relevant parameter to consider is R 3 uj 2 /(2GM p ) = 5 <C 1. One can then express <jW. as 



GM p 1 



»cff. = -- 



t^V +y) (27) 



2 

GM p GM P /~ 2 ~2 T 1 „ , 2 , 2 



,GMp GMp 
' R R 2 



,2 ,2 

X ■ V +z> 2 



(1 + 5)2 (1 + (5)2 



1/2 



which is exact up to terms of first order in -q and 5. Clearly, (f> c g is a function of r alone, where r 2 = 
J2i x i l a i 2 i where a\ = a 2 = 1 + 5, a 3 = 1. In the above derivation a\ = a 2 ^ as. A rotating body of 
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self-gravitating fluid can actually be triaxial in general, but triaxiality is likely unimportant for the rotational 
velocities of interest here (see Chandrashekhar 1969, Bertotti & Fasinella 1990). 

Putting the above into equation (p5[), we obtain 

p = p exp[-(r-Ro)/H] (28) 
where H" 1 is the derivative of the effective potential, and po is the atmospheric density at r = Rq. The 



choice of Rq is arbitrary at this point. In §2.1, we choose it to be the radius below which the atmosphere is 
completely opaque. 

Finally, note that the principle axes defined by x\ are not necessarily lined up with the axes defined by 



Xi in §2.1 by the lensing geometry. The two sets are generally related by a rotation. Since there is freedom 
in rotating x\ and X2 (i.e. the lensing geometry only picks out X3 or the z direction), we can without loss of 
generality relate the two sets of coordinates by two rotation matrices: x' = R z ' (7) ■ R x < (/3) • x, where R z ' (7) 
is a rotation about x' 3 by angle 7, and R x '(/3) is a rotation about x[ by angle (3. Using a\ = 03(1 + ca) and 
a 2 = 13(1 + £b), and assuming ea and e B are small (in our above derivation, €a = £ B = S i.e. an oblate 
spheroid; our derivation below continues to work even if this were violated), we can write 

r 2 = z 2 + (1 - 2eACOs 2 7 - 2e B sin 2 7)a: 2 + 4(e A - e B )xy + (1 - 2eAsin 2 7COS 2 /3 - 2e s cos 2 7Cos 2 /% 2 (29) 

where z = [1 + 0{e Al e B )][z + 0{e A ,e B )x + 0(e A ,e B )y}. 

Now, recall that the quantity that we are interested in, the deflection angle 9d, is given by equation (||): 
9d = Vj_ ap{r)dz. One can clearly change the variable of integration from z to z: apo exp[— (r — 
Ro) I H]dz = apo exp[— (r — Rq) / H]dz where we have absorbed the slight change in multiplicative factor 
into a redefinition of po. Moreover, since z is a dummy integration variable, we could as well rename z — > z 
in equation (E£ 



Consider next the terms involving x and y in equation (29). We could easily perform a rotation to put 
them in the simple form: (1 — e)x 2 + (1 + e)y 2 if we absorb multiplicative factors into a redefinition of H 
(such multiplicative factors would affect z also, but they can once again be absorbed into redefinition of po). 
To be more specific, let us consider the important case where 6a — £ B - Then, the x and y terms reduce to 

(1 - 2e A )x 2 + (1 - 2e j4 cos 2 /3) 2 / 2 = [1 - e A (l + cos 2 /?)] [(1 - e)x 2 + (1 + e)y 2 ] (30) 

where 

e = e A {l - cos 2 /3) (31) 

For a random distribution of angle (3, we expect on the average e = It is unclear, however, if f3 should 

be randomly distributed. The solar system planets do seem to have rotational axes pointing in all kinds of 
directions with respect to their orbital planes. The quantity €a describes directly the shape of the planet: 
contours of constant density obey z' 2 + (x' 2 +y' 2 )/(l + ca) 2 = constant. In other words, oblateness, defined 
by the ratio (ai — 0,3) /{ax + 03), is given by e^/2 to the lowest order. One can view e as a kind of projected 
ellipticity or oblateness. 

To summarize, with suitable rescaling of z, H and p$, we have 

ap exp[-(r-R Q )/H}dz , r 2 = (1 - e)b 2 + (1 + e)b 2 + z 2 (32) 



— OO 



where we have equated x and y with impact parameters 61 and 62 as defined in §2.1 
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Appendix B 

We discuss here how to solve equation (|B|) and how to find the caustic and critical curve. Equation (||) 
has a form that is exactly analogous to elliptic potentials sometimes used in modeling gravitational lenses. 
The main trick for solving this type of lensing equation is taken directly from Schneider, Falco & Ehlers 
(1992), but we provide additional comments here for cases that require special attention. 

Equation (||) is a set of two equations for 9\ and 9 2 , given the source position 8g- Eliminating ip(u) 
from the two equations, one obtains 

o2 _ (1 - e)gjg| r „„v 

Ul ~ (1 + ej¥ s - 2e9] [66) 

This gives us 9 2 as a function of 9j, and substituting into the first component of the lensing equation leaves 
us with a single equation for 9): 

91 = 9} + (I- e)9)^{u) (34) 

where u is a function of 9\ and 9j (now function of 9j) as well. The lensing problem is therefore no more 
complicated than one with spherical symmetry. Once 9\ is solved for, 9] can be obtained from equation 



For the case of 9% = 0, the above solution always puts Of = 0. There is, however, another possibility 
when multiple images are allowed. Examining the analog of equation ( [34] ) for the second component, one 
can see that the other option for a vanishing 9 2 s is to have 

1 + (1 + e)ip(u) = (35) 

This gives a single equation for 9] if one uses the fact that 9\ = (1 + e)9g(2e)~ 1 . This is obtained from 
eliminating ip(u) from the two components of equation (ph, but noting 9 2 s = and assuming 9 2 ^ 0. 

There is another special case: 9g = 0. In this case, equation ([54]) tells us that there are two possibilities: 
1 + (1 — e)tp(u) = or 9\ = 0. The former is handled automatically with the computational procedure above. 
The latter requires more care: one can solve 

9 2 s = 9j + (l + e)9ji>(u) (36) 

for 9] by setting 9\ = 0. Doing so ensures all possible images are uncovered. 
To find the caustic and critical curve, we solve A^ 1 = using equation (|^): 

v = ^1 [1 + 2<ip + (1 - e 2 ) + (1 - c 2 )u 2 ^ + u 2 ^} (37) 

This gives us v for a given u. However, from the definition of v and u in equation (|^), it is clear only \v\ < u 2 
is physical. Therefore, the caustic or critical curve corresponds to solutions of v given u that satisfies \v\ < u 2 . 
Strong lensing occurs when such solutions exist. For each pair of v and u, it is easy to solve for 9\ from their 
definitions: 



Given the critical curve defined by 0}, one can then solve for the caustic using the lens mapping (equation (Bb). 
We note that the above two expressions were interchanged by mistake in Schneider et al. (1992). 
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Gas comp. 


A 1 x 10 5 cm 3 


Bi x 10 n cm- 2 


Hydrogen 


13.6 


7.7 


Oxygen 


26.63 


5.07 


Nitrogen 


29.19 


7.7 


Air 


28.79 


5.67 


Methane 


42.6 


14.41 



Table 1: Refractive coefficients for different gas compositions. Refractive index is n = 1 + ap, with 
ap = A\{1 + Bi/X 2 ). Values are given for gases at 15° C and 1 atm. unit of pressure. 
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Fig. 1. — A schematic ray-tracing diagram: the solid line joining O (observer), L (lens) and S (source) 
represents the light ray. L represents the point of closest approach in the transiting planet's atmosphere. S 
represents a point on the surface of the star. The deflection is exaggerated and the diagram is not to scale. 
We use the convention that 6d has an opposite sign to Oj and Os- 
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Fig. 2. — Geometrical set up of a planetary transit. The circular solid line represents the boundary of the 
star, while the ellipsoidal dotted line represents that of the planet. The solid square is the center of the star, 
which moves along the dashed line with respect to the planet. The horizontal dotted line defines the major 
axis of the planet. 7tr.imp. is the angle between the horizontal and dashed line, while #tr.imp. is the distance 
in angular units of closest approach between the planet center and the stellar center. 
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Fig. 3. — The solid and dashed lines delineate the separation of the strong lensing regime (existence of 
caustics) from the weak lensing regime, for the combination of parameters B and Rq/H . The solid line is for 
projected oblateness/ellipticity e = 0.1 and the dashed line is for a sphere, e = 0. The parameter B is defined 
as B — 2a(po/H)DhsDoh/(Doh + Dls) (equation (^|)). The quantity Rq/H is the ratio of planet radius 
to atmospheric scale height. We show values of B and Rq/H for HD 209458b (HD), where Z?ls Adl, 
and for Mars (M) and Jupiter (J), where Dls ^ £>ol- In all cases, observations at optical wavebands are 
assumed. We show with the symbol J' the parameters for Jupiter if it were observed at wavebands with 600 
times stronger absorption (e.g. in the UV, see text for details). The solid square denotes the fiducial model 
we study in this paper. The triangle, open square and cross are other models we also discuss. 
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Fig. 4. — Upper Panel: The dashed line (overlaps with solid line in the upper portion) shows the magnifica- 
tion as a function of source position, for source positions that lie on the major axis of the planet (the x-axis 
in the lower panel.) The dotted line shows the occultation kernel W as defined in equation (||). The solid 
line represents the light curve for an occultation by a solar system planet of a distant star, which happens 
to be moving along the planet's semimajor axis. See text for details. Lower Panel: The solid line shows 
the caustic, source positions where the magnification diverges, while the dotted line shows the critical curve, 
which is the image of the caustic. The four solid squares are images of a point source denoted by an open 
square situated just inside the caustic. The axis coordinates arc distances (in km) in the lens plane, i.e. they 
correspond to physical distances from the planet center. 
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Fig. 5. — A transit light curve for our fiducial model (denoted by a solid square in Figure [| with atmospheric 
lensing (solid curve) and without (dashed curve). F(t) is the normalized stellar flux (equation (|t])). The star 
here has i?* = 0.6i?©. 
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Fig. 6. — The solid line shows fit), the fractional difference between the flux deficit from a stellar light 
curve with and without planetary atmospheric lensing (equation (|l4|)). This is for the fiducial model denoted 
by the solid square in Figure a. The dotted line shows f(t) as defined in equation (^q), the fractional 
difference between the light deficit from a lensed fiducial model and an unlensed model with parameters of 
the unlensed model tuned to minimize this fractional difference: Rq = 35400 km, H = 300 km, i?* = 417600 
km, u* = 0.832, v* = -0.2475, p = 1.68 x 10~ 5 , e = 0.05, w = 2.9741 x l(T 13 s-\ 7t r .im P . = 49.5°, and 
^tr.imp. = 69696 km/^OL- 
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Fig. 7. — f(t) (equation (M)) for three different planetary models corresponding to the open square (solid 
line), cross (dotted line), and triangle (dashed line) in Figure ||. Parameters for the dotted line resemble 
those for HD 209458b. 
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Fig. 8. — dF/dt for the fiducial model with e = 0.05, denoted by the solid square in Figure ||. The lower 
panel is a magnified version of the upper one. The solid line is with lensing, and the dotted line is without. 
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Fig. 9. — The two curves show f e (equation (|i6[)), the fractional difference between the light deficit from 
a transit of an elliptical planet and a spherical planet, for the fiducial model denoted by the solid square in 
Figure |^. The solid line is for a transit with lensing, and dotted line for a transit without. 
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Fig. 10. — f e (equation (|l7|)), the fractional difference between the light deficit from an elliptical fiducial 
model and a spherical model with parameters of the spherical model tuned to minimize this fractional 
difference. The solid line is for a model with lensing, and the dotted line is for a model without lensing. The 
derived parameters of the spherical planet are Rq — 35230 km, = 0.8048, = —0.227 for the lensed case 
(solid line), and i?o = 35270 km and the same u* and for the unlensed case (dotted line). The dotted line 
is essentially the fractional difference in the light deficit of an elliptical compared to a spherical planet of a 
similar total area. See text for details. 
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Fig. 11. — /asym as defined in equation ( |l8| ) with (lower panel) and without (upper panel) lensing. The solid 
line is for e = 0.1, the dashed line is for e = 0.05, both for 7t r .imp. = 45°. The dotted line is for the same e's 
but with Ttr.imp. =0°. In all cases, the rest of the parameters are those of the fiducial model (solid square in 
Figure [j]). If the curves in the upper and lower panels are plotted together, they almost completely overlap 
i.e. light curves with and without lensing have an almost identical degree of oblateness-induced asymmetry. 
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Fig. 12. — The upper solid and dotted lines show f a (t) (equation (pl|)) for cases with and without lensing 
respectively. The lower solid and dotted lines show the corresponding f a (t) (equation (p^)) which is a 
minimized function. The small values of the lower two lines demonstrates that an appropriate Rq can always 
be chosen so that the step function model approximates a more realistic absorption model to high accuracy. 
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Fig. 13. — The effect of a planet crossing an Earth-sized star spot on a stellar disk. Top panel: An Earth-like 
planet transit. Bottom panel: A CEGP-like planet transit. 



